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I . INTRODUCTION 

Free Electron Lasers (FEL) may have an impact on the 
United States Navy, not just as a weapon, but in many fields 
of science such as Medicine and Physics. The tunability and 
high power output of a FEL makes it useful for a broad range 
of applications. 

As a weapon, the tunability allows the user to select 
the frequency that will propagate most efficiently through 
the atmosphere. Also, the FEL has the potential to produce 
the power level needed for use as a weapon. 

The tunability allows a doctor to select not only the 
wavelength, but also the pulse-repetition rate of the laser 
during surgery. The FEL is able to take advantage of 
resonant modes in a molecule to allow the doctor to 
precisely remove millimeter layers of skin, fat, muscle, 
cornea, or tumor with no collateral damage to surrounding 
tissue. The FEL "scalpel" may revolutionize the operating 
room . 

The FEL is also the most viable path for achieving a 
coherent x-ray source. There are several projects working 
towards a high gain FEL that can achieve 1 angstrom x-rays . 
The uses for coherent x-ray sources are numerous. One use, 
which has a direct effect on the Navy, is to probe the 
reactions that occur in the initial moment of a nuclear 
blast . 
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There are numerous other uses for the FEL. Almost every 
field of science is finding new applications for the FEL, 
from probing light interactions with biological functions to 
studying Quantum Electro-dynamics. 

This thesis covers several different FELs which have a 
wide range of applications. The second chapter of this 
thesis covers the background theory that will be used in the 
rest of the thesis [reference (1)] . In Chapter III I 
calculate analytically the height of a tapered separatrix. 
In Chapter IV I analyze various high gain FEL's and use some 
dimensionless parameters to develop a simple model for each. 
The final chapter I simulate the Darmstadt FEL using a 
computer model and analyse the results. This analysis 
consists of simulations at various desynchronisms that 
determine power and gain expected in the FEL. There are 
instabilities that may be experimentally observed as well. 
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II . BACKGROUND THEORY 



The physics behind the Free Electron Laser (FEL) is 
based upon both spontaneous and stimulated emission from an 
electron. The two major components of the FEL are the 
relativistic electron beam and a static alternating magnetic 
field. The relativistic beam of electrons is accelerated in 
the transverse direction by the magnetic field. This 
acceleration initially causes spontaneous radiation from the 
electrons. Due to relativistic effects, the emitted light is 
radiated into a cone along the axis of the laser. The key to 
an operational FEL is the interaction between the electrons 
and the light field. The kinetic energy of the electron beam 
is converted to the radiation energy of the optical field. 
The electron pendulum equation and optical field equation 
describes the coupling. 

A. INITIAL CONDITIONS 

A classical approach has proven to be the most fruitful 
and experimental results agree with the classical 
results . [ref erence 1] 

1. Undulator Field 

The magnetic field of the FEL is created by a chain of 
fixed magnets, called an undulator. There are two basic 
types of undulators : helical and linear. The linear 
undulator is a series of bar magnets with the polarization 
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in alternating directions, so that an alternating magnetic 
field is set up. Figure 1 depicts a linear undulator with 
the z direction is along the axis of the laser. The magnetic 
field for the linear undulator is given by 

B = [b^, bJ = bJo, sin(kcz)cosh(koy), cos(koz)sinh(koy)], (II-l) 

where k^ = 2n / is the undulator wavenumber. 



mirror 




Figure II-l Linear Undulator (after reference 7) 

The helical undulator consists of two concentric 

currents that create a magnetic field that rotates spatially 
along the length of the undulator. For the helical undulator 
the magnetic field is given by 

B = Bjcos(koz),sin(koz),0] (rj.o) 
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2 . Initial Optical Field 

To begin analyzing FEL operation we need to assume that 
the FEL is already operating and there is already a 
radiation field present. The optical field for the linear 
undulator is assumed to be a plainer, plane parabolic wave 
with the electric and magnetic fields given by 

Es = e[cos i|/,0,0] and Bs = e[ 0, cos n/,0] respectively, where 
v|; = kz - cot + (j) . The electric field amplitude E is given in 
cgs units with ^ being the optical phase and k = 2n / X is 

the optical wave number, co = 27 tc / X is the optical angular 
frequency . 

The initial optical field for the helical undulator is 
assumed to be a circularly polarized plane-wave. Hence the 
electric and magnetic fields are given by 

Es = e[cos v|;,- sin v|/,o] and Bs = E[sin \|/, cos i|/,0] , respectively. 

B . ELECTRON DYNAMICS 

1. The Electron Phase 

Before moving on, we need to define the electron phase. 
The electron phase is the key to coupling the optical field 
with the bunching of the electron beam. In the combined 
undulator and light field the electron phase is defined as 

^ (t) = (k+k(,) z (t ) +(0t . The electron phase definition involves 
two fixed wave numbers k = 2% / X , describing the light wave 
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and ko = 2n / describing the undulator, with the only 

variable being z(t) . The electron phase follows the electron 
position. The time spent by the electron in the undulator is 

At=L/(3()C»L/c~10'®s , where L is the length of the undulator 

and the speed of the electron (V 2 =cp 2 s»c) is approximately the 
speed of light. This time is extremely small, so it is 
easier to introduce a dimensionless time x and a 
dimensionless time derivative 



The dimensionless time x ranges from 0 at the beginning of 

the undulator to 1 at the end of the undulator (0<x<l) . The 
electron phase velocity is then defined as the derivative of 
the electron phase 



To study the operation of the FEL, it is most 
informative to plot the phase velocity versus the electron 
phase in what is called phase-space. This will be discussed 
further in subsequent parts of this chapter. 

2 . Wavelength of Light in FEL 

The traditional laser requires a bound electron to 
transition to a lower energy level and radiate. In the FEL, 
a relativistic beam of free electrons amplifies the light 




(11-3) 



V (X) = ^ = ; = L[(k + ko)p, - k] 
dx 



(II-4) 
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wave. The undulator wavelength A.^ is Lorentz contracted to 

r 

^0 = ^0 / Y electron frame. The emitted optical field 

is Doppler shifted is also 

V = (l - p,)yX « 2yX . (II-5) 

At resonance, there is optimal energy exchange. When 

t 

Xg = X' in the beam frame. When converted to., the .lab, frame 
the resonance condition becomes 
. . 1 + 

X = Xg ^ (II-6) 

2y 

where the dimensionless undulator parameter K is defined as 



K = 



eBA, 



eB 

kgmc" 



(II-7) 



2tx me 

B is the RMS magnetic field, e is the magnitude of the 
electron charge, m is the mass of the electron and c is the 
speed of light. 

3 . Electron Paths 

To analyze the electron trajectories in the PEL, we 
will look at the linear undulator. The magnets are aligned 
along the y direction. An electron that moves off axis in 
the y-direction feels a stronger magnetic field because it 
moves closer to a magnet . The electron is pushed back toward 
the axis creating a natural focusing. The focusing of the 
undulator causes the electrons to wiggle over many undulator 
periods in the focused direction. The wiggle is called a 
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betatron oscillation. To analyze betatron effects, we will 
consider a perfectly injected electron in a linear undulator 
that focuses in the y direction. Perfect injection means 
that the electron is injected into a perfect sinusoidal 
trajectory with no angular spread. 

The three spatial components are of the Lorentz force 
equation with the magnetic field for linear undulator from 
equation (II-l) as follows 

— [p, sin(koz)cosh(koy) - p cos(koz) sinh(koy)] , (II-8) 

ymc 

p = — [Px cos(koz) sinh(koy)], (II-9) 

ymc 

p, = - — [Px sin(koz) cosh(koy)] . (II-IO) 

ymc 

In the next section, we will derive these equations in 
greater detail for the helical undulator. The P^ term has an 
exact solution, and perfect injection implies the constants 
of integration are set equal to zero, so 



P. 



V2K 

Y 



cos(k(,z) cosh(kpy) . 



(II-ll) 



Substituting P^^ into the other two equations leads to 
r-K^k 

P = cos-(koz) sinh(2koy), 

r 



( 11 - 12 ) 



and 



P. 



cK^k„ 



sin(2koz) cosh’(koy) . 



( 11 - 13 ) 
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Since the y and z components are coupled, we add the 
restriction of perfect injection in the y direction and 
analyze the motion in the xz-plan. Hence, 

Px = cos(koz) (11-14) 

y 

py(0) = y(0) = 0 (11-15) 

cK^k. . . X 

= 7-2- sin(2koz) (11-16) 

r 

In the typical operating PEL K«1 and y>>l, so that 
K/y<<l. Expand in powers of (K/y) , so the x and z components 
become 



/ X 

z(t) = cp,t - - — 7 sin(2k(,cp,t) -i- . . . , 



4k(,y 



(11-17) 



/ X 

x(t) = sin(koCp^t) -t- . . . , (11-18) 

koY 

where P^ is the average velocity of the electron beam in the 
undulator. Then converting to dimensionless time x gives 



z(x) = p^Lx 



^ sin(47:Nx) -i- . . . , 
87ry 



(11-19) 



x(x) = — sin(27iNx) -t- . . . 
V27ty 

From these two equations, 
the transverse oscillation 



( 11 - 20 ) 

we can see that the amplitude of 
is / V27ty where the electron 
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oscillates once every undulator period. In the z direction, 
the electron oscillates twice every undulator period. 

Having looked at the oscillations in the x and z directions, 
we now go back and assume non-perfect injection to see what 
happens in the y direction. The oscillations in the y 
direction are called betatron oscillations. Assuming small 
oscillations so that koy<<l and sinh(2kpy) w 2k,^y we.,, can 
rewrite equation (11-12) and convert to dimensionless 
parameters to get 



Averaging over several undulator periods, since one betatron 
oscillation is typically over many undulator periods, we 
have 



Equation (11-22) is the equation of motion for a simple 
harmonic oscillator. The solution is a well-studied 
sinusoidal wave, with the dimensionless betatron frequency 

C03=KLk;,/Y=27tNK/Y. The number of betatron oscillations along 
the undulator is 



oo ^ 

y(x) = ^ cos^(27tNT)y(x) 



y 



( 11 - 21 ) 




( 11 - 22 ) 



NK 



(11-23) 



N 



Y 
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A typical value is Npwl, since N«100, K=l, and y=100 are 

typical values. FELs with larger values will be discussed 
later . 



C. ELECTRON PENDULUM EQUATION 

The dynamics of the FEL may be described very simply by 
the electron pendulum equation. The first step in deriving 
this equation is to consider an operating FEL where an 
optical field is already present. This means the electron 
experience forces from both the optical field and the 
undulator field. 

The helical field is much simpler to work with 
mathematically and the final formulas can be converted 
easily to the linear case. Hence, we will just use the 
helical magnetic field. 

The forces from the undulator field are governed by the 
Lorentz force equations (cgs units) , 




(11-24) 



and the energy equation 




= -e cP • E , 



(11-25) 



dt 



where 




(11-26) 
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The momentum is p=Ymv=ymPc, its energy is ymc^, and the 



charge magnitude e=|e|. We have five equations and four 



unknowns x, y, z and y. 

Substituting the magnetic field from equation (II-2) 
into equations (11-24), (11-25) and (11-26) yields the 
following equations 

e 



dt 



dt 



me 



[e( 1 - P,Xcos sin vj/,0) + 
P^b(- sin(koz), cos (koz),0)] 



= — — [e(Pj. cos v|/ - Py sin v|/) + 



me 



b(Px sin(koz) - py cos(koz))] 



d(y) 

dt 



= - — e[p^^ cos - py sin 1|/] , 



me 



(11-27: 



(11-28) 



(11-29) 



where 



P_L = P^ + Py and Vj/ = kz - cot + (j) . There are more 



equations than unknowns, so we can ignore the P^ component. 

The electrons are traveling near the speed of light 
along the z axis, so we can estimate that P^~l . This means 

that e(i - PJ « BP_ , and the first term of the transverse 
components of (11-22) can be ignored. Thus, 



s - — [P,b(- sin(kpz), cos(koz),0)] , 
dt me 



(11-30) 



Which implies that the static undulator field controls the 
transverse motion of the electron almost entirely. We 
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can integrate (11-30) by inspection, yielding the following 
results 



— K 

Pa = [cos{kgz), sin(koz),0] . 

y 

Substituting (11-31) into (11-29) yields 
eKE 



(11-31) 



y = 



ymc 



cos(^ + (()) 



(11-32) 



where ((^ + (|) = \|/ + k^z = (k + kQ)z - cot + cj) . From (11-31) 



2 K' 



P/ = — . 

y 

and substituting into to (11-26) 

1- = 1 - (p)* = 1 - - p/ = 1 - p.= - 



(11-33) 



(11-34) 



so that 



(i + = 1 - P3' • 

y 

Taking the time derivative of (11-35) , 

i = y'PzPz 
y 1 + 

From the derivative of the electron phase, we have 

4 = (k + kg)z - (0 = (k + kQ)cP - CO , 

and 

i; = (k + ko)cp, , 
so that 



(11-35) 



(11-36) 



(11-37) 



(11-38) 



13 



(11-39) 



Pz = 



c 



(k + ko)c ■ 

Substituting into (11-36), we get 
i ^ Y^PzPz ^ 

Y 1 + (l + K^Xk + ko)c ■ 

Now equating (11-32) and (11-40), we get 



Y^PzC 



eE,K 



(l + K^)(k + ko)c me 



cos (i; + (j)) , 



and solving for , we have 
_ (l + K^) (k + ko)c eEpK 



; = 



cos ((; + (!)). 

Pj Y me 



(ii-4o: 



(11-41) 



( 11 - 42 : 



since k>>ko_ we can use the approximation (k+k(,)«k. From 

equation (II-6) cOg = ckg = ck(l + K^)/ 2y~ , and using 

equation (11-42) reduces to 
•• 2c0geEK 

C = — I cos(C + (|)) (11-43) 

Y me 

Converting to dimensionless time with dt=Ldx/c, and defining 
the dimensionless optical field magnitude as 



a = 



4:i:NeKLE 

2 2 
Y me 



The pendulum equation becomes 



(11-44) 



= |a| cos ((^ + <t>) . (11-45) 
The pendulum equation (11-45) is the same form for a linear 
undulator. However, for the linear undulator, the optical 
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field |a| is modified by substituting K [ (^) - (^) ] for K, 
where ^=K^/2(1+K^) and Jq (^) and Ji (^) are Bessel functions. 



D. OPTICAL WAVE EQUATION 

The pendulum equation describes how individual 
electrons are affected by the laser field. To describe how 
the laser field is affected by the beam of ,, electrons, we 
next derive the optical wave equation. Begin with the 
inhomogeneous wave equation for the optical vector potential 

A, 







(11-46) 



where 



E 



S 



1 ^ 

c dt 



e[cos \\i,~ sin , 



(11-47) 



Bs = V X A = [sin v(/, cos v|/,0] , (11-48) 

i|y = Icz - (Ot + <j) . When the amplitude and phase of the laser 
beam change slowly with time, we can write the vector 
potential in the form 

A(z, t) = — E(z, t) [cos (v|/) ,- sin (y) ,0] . (11-49) 

CO 

The form in equation (11-49) assumes that the electron beam 
is small compared to the waist of the optical mode, so that 

the transverse direction is not important . A is dependent 
only on z and t. The spatial derivative is 
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a"A 1 dE 



dz~ 



k dz 
1 d^E 
k dz^ 
1 ^ 
k dz 

52 



5(1) r 1 

k H [cos (\|/) sin (\|/) ,0J + 

dzj 

[sin (v|;) , cos (i|;) ,0] + 

[cos (v|y) sin (v|/) ,0] + 



5f 
k + — 
dz 



ill-50) 



E 5 d) r T 

- — T [cos (v|;) sin {\\i) ,0J 
k dz 



+ 



E 

k 



5z 



- 2 

[- sin (if/) cos (If/) ,0] 



and the time derivative is 
CO d^h dE r 5(|) 



c 5t’ 5t 



5t 



- CO 



5^E 



— 

~dt 



5(1) 



- CO 



[cos (i|/) sin (i|/) ,0] + — r [sin (i|/) , cos (if/) ,0] 

5t 

5(1)' 



[cos (i|/) sin (if/) ,0] + E — r [cos (if/) sin (i|/) ,0] (II -51) 

5t 



+ E 



5t 



- CO 



[- sin (ij/) cos (if/) ,0] . 



From the assumption that phase and amplitude varies 
slowly, 5 e/5z << kE , d^/dz << k(l) , 5E/5t << coE , and d^fdt. 

<<co(l) . From these approximations, we find that 



1 5' ' 




'5E ^ 1 5E‘ 


c' 5t' 


A ~ 2 


_ 5z c 5t 



[cos(i|;),sin(i|/),0] 



(11-52) 



+ 2E 



d^ 1 5 ( 1 ) 

5z c 5t 



[-sin(i|/) -cos(i|/),0] = -— . 

c 



To further simplify, we use a coordinate transformation, 
where the position z moves along with the optical field at 
the speed of light c, so that z = z + ct . Also, we convert 
to dimensionless time, t=ct/L to get 
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(11-53) 



A_1 - A_ - A_ 

ydz cdt) cdt Ldx 
and 



lA. 

L dr 



[cos(v|;), sin(y),0] + 2E 



L dr 



[-sin(v|/),-cos(v|/),0]^ 




(11-54) 



Each electron adds to the total current, so the current 
density contribution from the i*^^ electron is summed to get 

=-ec^p^5^^^(x-rJ (11-55) 

i 

where the Tj is the trajectory of the i"^*^ electron and 



is a three dimensional Dirac Delta Function. Substituting, 
the transverse velocity = -K / y(cos(k(jz), sin(lcQz),0) into the 
current gives 

- — ~ ——[cos (k(jZ) , sin (k^z)] -ri). (11-56) 

c i Y 

Therefore, the wave equation is 



5e' 

dr 



[cos (i);) , sin (i)/) ,0] + E 



dr 



[- sin (V|/) ,- cos (V|/) ,0] = 



^ _ 27ieKL ^^^^ ^ g(3)|^ - rj 



(11-57) 



Both sides of (11-57) averaged over a volume element dV^, 
which is smaller than the coherence volume, yet larger than 
a wavelength of the light. The current density over this 
small volume element is constant, so the summation over i 

electrons ^ (- 2neKL / y)[cos (kj,z) , sin(kQZ)] can be replaced 
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with a charge density p times the average of 

electrons in the beam, the average is symbolized by 

now we have two uncoupled equations 

dE 27repLK , . 

— = — (cos (C + <t>)) 

dr Y 

d(j) 27iepLK , . > 

— (sin(i; + (j))} 

ox y 



sample 
{ ) . So 



(11-58) 



(11-59) 



where i|;+koZ=(I^+(t). Defining the dimensionless current 
j = 8t'e~pNK^L^/y^m c^ , and recalling the dimensionless optical 

field |a| = 4TcNeKLE/Y‘itc , the change in electric field 

amplitude becomes 

dial “ > 

= a = -j(cos (C + (11-60) 

dr 



Since the electric field is complex, then a = |a|e^* , and we 
arrive at the optical wave equation 



a = -j(e"^). (11-61) 

Equation (11-61) is derived for a helical undulator. 
However the same form of the equation applies to a linear 

undulator provided K is replaced by K [ Jq (^) - (^) ] . 



E. GAIN MECHANISMS 

For the FEL to be operational, the there must be an 
increase in power of the optical field each pass through the 
undulator. There are two regions of optical gain: the low 
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gain region {j<7r) and the high gain region (j>> 7 t). To 

understand the mechanisms behind gain we first develop the 
concept of an electron phase space plot. 

1 . Phase Space 

The electron phase space is a plot of an electron's 

O 

phase velocity, v(i) = (^ = L[(k + - k] and the electron 

o 

phase (t) = (k+ko) z (t) +cot . In the low gain case a w 0 , which 

means that equation (11-61) evolves slowly so that |a| and ^ 
are approximately constant. Multiplying both sides of the 

O 

electron pendulum equation (11-45) by yields 

O 0 o o o 

c = |a| cos (^ + (})) (^ = V V (11-62) 

/ 

or 



V V - |a| cos ((^ + (|)) (^ = 0 
so that 



(11-63) 



dx 



y - |a| sin((; + (f)) 



= 0 



(11-64) 



This means that - 2|a| sin((!^ + <})) = cons tan t , and 

V- = Vq- - 2|a|[sin((; + (t>) - sin(Co)]. (11-65) 

Where Vp=v(0) and ^o=(^(0) are the initial phase velocity and 
phases. Equation (11-65) is a statement of energy 
conservation in the dimensionless coordinates and v. The 
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equation also identifies paths in the phase space, v {Q , 
that describes the electron motion. In Figure II-2, several 
phase space paths are plotted. 5 refers to taper and will be 
discussed in the next section. 



4Va 



V 



Figure II-2 Phase Space Diagram 
2 . Separatrix 

Electrons stay on v(d^) paths in the plane. These 

paths are based upon the initial phase velocity Vg and 

initial phase of the electron and the pendulum equation. 
This same equation describes an undamped, undriven pendulum 
where is the tilt angle and v is the angular velocity. The 





PHASE SPACE 


PATHS. . 


a=l 


(j)=0 


5=0 
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critical points for these paths occur at v=0 and t^=%/2±nn 

where n=l,2,3. The stable fixed points occur when C^=% / 2±2nn . 
In Figure II-2 the stable fixed point is labeled with an 
"S" . The unstable fixed points are at v=0 and d^=37i/2±27in . In 

Figure II -2, the unstable fixed points are labeled with a 
"U" . At these fixed points, the electron remains stationary- 

in phase space. At {v, (^) = (ti/2 , 0 ) , the pendulum is hanging at 

the bottom of the arc. When iv|>2|ah^^, the pendulum has 
enough angular momentum to swing continually over the top 
around in a circle. In phase space, we can draw a separatrix 
which separates the lower closed orbits from the over-the- 
top open orbits. The equation for the separatrix is 

= 2 |a|[l-sin((^ 5 +([))] . (11-66) 

In Figure II -2, the separatrix separates the open orbits, 
shown as a gray area, from the closed orbits, shown as a 
white area. 

3 . B^mching 

Consider a collection of electrons that are initially 
uniformly distributed in phase evolving along the undulator. 

OO 

The pendulum equation shows that for -7t/2 <(,<%/ 2, <^>0 and 

the electron phase accelerates, while for 7i/2 < < 3n/2 , 

00 

(^<0 and the phase decelerates. The result is that the 
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electrons tend to bunch about 7i/2 as they evolve along the 
undulator . 

4. Changes in Phase Velocity 

From the definition of phase velocity 

v(t) = L[(k + k„)p. - k], we find the change in beam energy due 
to a change in phase velocity. Taking the derivative of both 
sides near resonance where v=0, k>>ko and N=L/A,o, we get 
6v w LkSf^z NX^kSp^ . From equation (II-6) and the resonance 
condition, A = Aq(i + K^)/2y^ / we find that 
c 5y 

8v == 4 kN — w 27 tN — . (11-67) 

Y 

Thus, an increase in v results in an increase of the 
electron energy, ymc^, and the optical wavelength X. 

5 • Gain Formula 

Gain is defined as AP/P, where AP is the change in optical 

power in one pass through the undulator, and P is the power 
at the beginning of the undulator. From (11-67) the average 
change in energy for an electron in the undulator is 

Aymc" w ymc^((v) - Vq)/47xN, (11-68) 

where (v) is the average final energy. The number of 
electrons in a volume dV of the optical wave is given by 
pFdV where p is the electron density and F is the filling 
factor. The filling factor F is the cross-sectional area of 
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the electron beam divided by the cross-sectional area of the 
light beam. The radiation energy in a volume dV is given by 



2E^dV / 8 tc in cgs units. Invoking the definition of the 
dimensionless field a and the dimensionless current j , we 
arrive at 



G 



[pFdv][yrac^((v) - v^) / 47 in] 

[2EMv / 87c] 




(11-69) 



6 . Low Gain 

For low gain, the optical field changes vary slowly so 



we can estimate that a w 0 with the initial conditions 
v(0)=Vq_ a(0)=ao, and C^{0)=C^^. From perturbation theory, we 
find that 



V = Vq + -% [sin(;o + v,x) - sin((;J] + 

Vo 

- \ cos(2^o + 2V(,t) - cos(2^o) + cos(2v,x) - 



(II-70: 



+ . . . 



_1 - v,x sin((;J cos((;o + Vqx) 

Taking take the time average of the phase velocity, we find 



(v) = Vq + -% [2 cos(vqx) - 2 + VqX sin(vox)] + . . . 

Vq 

The gain is therefore 
2*i F* "i F 

G = — (vq - v> = — [2 - 2 cos(vox) - VqX sin(vox)] . 



(11-71) 



(11-72) 
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7. 



High Gain 



In the high gain regime the current is very large, 
j>>m. Unlike the low gain regime, we cannot ignore the rapid 

evolution of the optical field over a single pass. Combining 
the expanded pendulum and optical equation and integrating 

over all the initial electron phases, Jd^o , we get 

a(x) = ^ JdTYF(x')e'^''°"'a(T - x') (II-73) 

where F(t')= Jdqf(q)e‘'‘'^ is the characteristic function of the 
distribution f (q) . f (q) is the distribution of the initial 
electron phases v^=Vg+q about Vg and |dqf(q) = l. 

The low current gain formula can be obtained from this 
equation by considering a perfect beam where f(q)=5(q) and 
a(x — x') « Sg , in which case a^ can be moved outside the 

integrand. Further integration results in the low gain 
formula shown earlier. 

To find the solution to equation (11-73), we need to 
consider a perfect beam with F(x')= 1 and a resonant optical 
field Vg=0. Then, taking two derivatives, we get 

a°(x) = ^a(x). (11-74) 

It can be shown that the solution to (11-74) is 
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(11-75) 




/2 ^ g(3/2)‘' = (i-V3)x/2 g( 






Since for higli gain j>>7t the dominate exponential terms are 
the fastest growing terms. The following approximations can 
be made for the optical field amplitude, phase and gain: 



The gain grows exponentially along t with a growth rate 



F . BEAM QUALITY 

Previously, when considering the interaction of the 
electron and optical pulse we assumed perfect injection. In 
reality there are errors involved. The two factors that have 
the most influence are electrons injected off axis at the 
initial position, yo, and electrons injected at an initial 

angle, 0y. From Reference (2) the decrease in phase velocity 
is 




(11-76) 




(11-77) 




(11-78) 



proportional to 



27tN 




(11-79) 
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The change in phase velocity from these errors is made more 
severe with more undulator periods N. A longer FEL is more 
sensitive to injection errors. 

A measure of the quality of the electron beam is the 

emittance, = YoSy , where is the rms initial position 

spread of electrons along y, and 0y is the rms initial 
angular spread of the electrons in the y direction. The beam 
can be focused to minimize either Yq or 0^ , however 

emittance remains fixed. The minimum change in phase 

velocity occurs when the beam is "matched" so that both 
terms in (11-79) contribute equally. Matching occurs when 
Kkoyo ~ y9y . For a matched beam Av « 27tNoKk.Sy / y . 

To show the effects of beam quality we need to look at 
the FEL integral equation (11-73). The function f (q) in 
equation (11-73) has a distribution based only on the 
initial distribution of the electrons as they enter the 
undulator. Two examples of distributions are those from the 
electron phase velocity and injection angle. The spread in 
phase velocity can be considered a Gaussian distribution 

about Vji with a standard deviation of cr^ = 47iNAy / y . The 
angular spread leads to a exponential distribution with a 
standard deviation given by Cg = 47tNy^0y^ /(l + K^). The 
distributions are 
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(11-80) 



and 



fe(q) = 



e 



(11-81) 



The characteristic function for these two cases are 




(11-82) 



and 



1 



Fe(q) = 



(11-83) 



O ' -i' . ' ' 

1 - IOqT 

For perfect injection ag=0 and ae=0, we have F(x)=l. When the 



decays with a characteristic time l/a^ or I/oq. F (t) depends 

on the shape of the function f(q), and different 
distributions of f (q) effect the FEL's ability to bunch. 

G. SHORT PULSE EFFECTS 

In the low gain FEL, the modes of light in the 
oscillator can develop independently. Since the electrons 
are relativistic, they travel to, but less that, the speed 
of light. At resonance, one wavelength of light passes over 
an electron as the electron travels through one undulator 
period. Over the total length of the undulator the electron 

slides behind the light by NX, called the slip distance. In 
the transverse dimension, the optical mode is roughly the 



beam is not perfect the characteristic function 
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radius of the optical beam which is (NX,/7 i) . These modes 

evolve as determined by the pendulum equation and wave 
equation . 

In order to follow the evolution of the optical field, 
we must set up multiple sites in z along the complex wave 

envelope. The transformation to spatial modes, a— >a(z), is 
the equivalent to the transformation to wave numbers, 
a->a(k) . The longitudinal position on the z-axis Z is 
transformed to dimensionless parameters by dividing by the 
slippage distance Nl, therefore z->Z/N>,. Througout the 
remainder of this paper z is dimensionless parameter. 

The length of the electron pulse Ig becomes the 

dimensionless parameter CT 2 =le/NA,. The pendulum and wave 

equations of motion with slippage included become 

oo 

<^z-T =K|C0S(C,_, +({)J (11-84) 

and 

■ (11-85) 

The electron current is = 87z^e"p(z)NK^L^/y^mc‘ , where the 

electron density, p(z), is the density at the site z. A 

site in the electron beam interacts with and amplifies N 
optical wavelengths, while the optical wave passes over and 
attempts to bunch electrons in a slip distance along the 
electron beam. 
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Short pulses of spontaneous radiation drive a typical 
FEL oscillator. When the length of the electron pulse, 1^, 
is comparable to the slip distance, short pulse effects 
influence the interaction. As the electron and optical 
pulses pass over each other, the pendulum equation (11-85) 
sees a continually changing optical field, a^, while the 
wave equation (11-86) sees a continually changing current 

density j^.t and bunching distribution . In the 

resonator, the optical pulse bounces between two mirrors. 
The distance between the mirrors is S which is necessarily 
greater than the length of the undulator, L. The time 
interval for the light-pulse make one round trip is 2S/c. 

Desynchronism, d, is the difference between the arrival 
times of the electron and optical pulses at the beginning of 

the pass, again normalized by the slip distance NX. If d=0 
at T=0, the electron pulse arrives simultaneously with the 
light-pulse. Moving one of the mirrors slightly by AS 
changes d so that d=-AS/NX. 

H . TAPER 

1 . Purpose of Taper 

In a successful FEL design, the optical field will 
eventually saturate due to over-bunching of the electrons 
along the undulator, thereby changing the dynamics of the 



29 



electron interaction. As the electron loses energy to the 
optical field, the phase velocity moves off resonance by 

Av » 4TrNAy/y « -2 ti . One way to extend saturation to higher 
powers is to taper the undulator field. Instead of a 
periodic sinusoidal magnetic field with constant amplitude, 
the amplitude of the field is either increased or decreased 
progressively along the -length of the undulator. Either 
shortening the undulator period along the undulator, or 
changing the strength of the magnetic field along the 
undulator tapers the undulator. In either case, acceleration 
of the phase of the electron results . 

This acceleration is denoted by 5, or the taper 

O 

parameter, so that v=5, and 

kak 

5 = -4tcN ^ , (11-86) 

1 + 

where the undulator field is tapered by AK over the 
undulator length. The electron pendulum equation becomes 

oo o 

(;^ = v = 6 + |a|cos(^ + (j)) . (II -87) 

While the optical wave equation (11-86) remains unchanged. 

2 . Taper Phase Space 

In the limit of low gain the phase space paths for the 
tapered undulator are given by 

1 

V = [vq^ - 26(C - Co) - 2|a|[sin(C + (j>) - sin(Co + <t))]] ^ (11-88) 
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The stable fixed point for both positive (6>0) and negative 
taper (5<0) in phase space is at 

(^, u) = (cos'' (- 6/|a|)- <(),o) (11-89) 

The unstable fixed points depends upon the direction of the 
taper. The unstable fixed point for 6>0 in phase -space is 
located at 

(C, u) = (271 - cos'' (- 6/|a|)- (}),o) (11-90) 

and for 6<0 it is at 

((;,u) = (- cos”' (- 6/|a|)- (|>,o) . (11-91) 

The formula for the separatrix is 

- Co ) + 2|a|[sin((;, +(}))- sin((;o +<{))], (11-92) 

where (^o=27t-cos'^ ( -6/a) for 6>0 and <!^()=-cos‘M -6/a) for 6<0 . 

Figure II -3 shows phase space plots for negative and postive 
tapers. "U" refers to the unstable fixed point, and "S" 
refers to the stable fixed point. 
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Figure II -3 Phase Space Plots for Negative and Positive 
Taper 
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Ill, TAPER SEPARATRIX HEIGHT 



To find the height of the separatrix for the tapered 
undulator it is convenient to make a few variable changes 
and assume that there is no initial phase (which only shifts 
the separatrix along x) (!^ = x, v = y, (|) = 0, and 

= Xq . Hence, equation (11-95) becomes 
y^ = 25(x - Xg) + 2|a|(sin x - sin x^) (III-l) 




Figure III-l Separatrix for negative taper with q=0,4 

Figure III-l shows the separatrix for a taper 5=-2.5 and an 
optical field ao=20. The height of interest is the maximum 
height labeled h in Figure III-l. To find the maximum value, 
we need to take the derivative of equation (III-l) which 
gives 

2yy' = 25 + 2|a| cos x , (III-2) 
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where y' = dy / dx , and set the slope y' equal to zero 
5 + |a| cos X 

y = = 0. (III-3) 

y 

The maximum value for x can be denoted as x and is given by 
5 + |a| cos X = 0 or x = cos'^(- q) , (III-4) 

where q = 5 /|a| . Note that q is restricted to values between 
1 and -1. 

A. POSITIVE TAPER 



For 5>0, Xo=2ti:-cos‘^ ( -q) and substituting x for x gives 

y^ = 25[cos"^(- q) - izn - cos”^(- q))] + 

r / \ / M (III-5) 

2|a|[sin(cos'^(- q)j - sin(27i - cos ^(- q))J 

Using the identity sin(2Tc - x) = -sin(x) we can simplify equation 
(III-5) to get 

y^ = 25(2 cos"^(- q) - 2n] + 2|a|[2 sin(cos'^(- q))] . (III-6) 

Factoring out |a| we have 

y^ = 4|a|[q(cos'^(- q) - m) + sin(cos'^(- q))] . (III-7) 

Using the identities cos“^(- x) - % = - cos’^(x) and 
sin(cos"^(- q))= VW. we can consolidate equation (III-7) 



to get 
y' = 4|a| 



- q cos 



(q) + Vl - • 



(III-8) 
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Taking the square root of both sides, we get the maximum 
value of y to be 



y 



= ^ cos"'(q) + 



(III-9) 

The height of the separatrix, h is twice the value of y 
therefore. 



h = ^Vn[" q cos'^(q) + Vl - 



(TIT -10) 



B . NEGATIVE TAPER 

For negative q, Xq = - cos”^(- q) . Therefore equation 
(III-l) becomes 

y^ = 25[cos"'(- q) - (- cos“^(- q))l + 

r / ^ / M (III-ll) 

2|a|[sin(cos ^(- q)j - sin(- cos ^(- q))J . 

Utilizing the identity sin(- x) = - sin(x) , it becomes 

y^ = 45[cos'^(- q)] + 2|a|[- sin(cos'^(- q)) - sin(cos"^(- q))] (III-12) 



Consolidating and factoring out |a| gives 

y^ = 4|a|[q(cos'^(- q)) + sin(cos'^(- q))] . 

Take the square root of both sides gives 



(III-13) 



- 1 



y = 2a/Rr cos'^(- q) + 



and then the height is 



h = 4^_q COS'H- q) + yjl - q^] 



(III-14) 



(III-15) 
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Equations (III-IO) and (III-15) are equivalent with the |q|. 
Therefore we can say for all values of q 



h = N cos'^(jq|) + 



(III-16; 



C. POSTIVE AND NEGATIVE TAPER 

Since the separtrix height for the non-tapered case 
(5=0) is 



h = " N cos'"(jq|) + 



(III-17: 



Figure III-2 shows a graph of equation (111-18) . 




Figure III -2 Change in Separatrix Height 5/a 
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D . LINEAR APPROXIMATION 

A Taylor expansion inside the brackets of equation 
(III-18) and using the expansion of arccosine 

cos‘^(q) = — - q + . . . (III-18) 

2 



gives 



h = 



n 

T - q 



- q 

which simplifies to 

1 ^ 1 2 

1 q + - q 

2 2 






h = 4VN 



(III-19) 



(III-20) 



h = 4VNI 1 - - q 



7t 






= 4Vn[ 



5^ 

1 - 0.785 - 

a) 



(III-21) 



From equation (III-21) it can be shown that, as 5 -> 0 the 
separatrix height h —> , or the size becomes the same as 

the untapered case. As the optical field a — > oo and the 
height h — > the same as 6=0. As the optical field 

becomes zero a — >• 0 then the height h ^ 0 . 
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IV. PARAMETERS FOR USE IN SIMPLE MODEL OF A HIGH GAIN FEL 



To develop a simple model for a high gain FEL we need 
to first take a close look at the parameters that dominate 
the electron interaction with the optical field in the high 
gain regime. Table IV- 1 is a table of several proposed and 
well -documented high gain FELs . 





Units 


LCLS 


TESLA 


ELF85 


ELF 8 6 


ELF 8 7 


MIT 


KE 


MeV 


1500 


35000 


3.3 


3.5 


3.5 


2.3 


I 


amps 


3400 


5000 


500 


850 


1000 


930 


X 


cm 


1.5x10'® 


1x10'® 


0.76 


0.72 


0.175 


0.125 


L 


m 


100 


60 


3 


3 


4 


1.5 


j 


- 


27200 


35000 


3300 


5000 


5600 


3800 




- 


0.45 


0.003 


0 . 045 


0 . 04 


0.12 


0.16 


^0 


- 


0 


0 


12 


14 


0.34 


0 




cm 


3 


5 


9.8 


9.8 


9.8 


3.14 


N 


- 


3300 


1200 


30 


30 


40 


48 


Np 


- 


0.3 


0.4 


11 


11 


5.7 


2.3 




- 


8.4 


3 . 7 


12 


12 


32 


3.6 


le 


cm 


0 . 008 


20 


450 


1498 


300 


900 


NX 


cm 


3x10'® 


2x10'® 


23 


22 


7 


6 


Slip 


- 


186 


1000 


19 


68 


42 


6,000 


Wox 


cm 


0 . 002 


0.006 


5 


5 


1.6 


0.8 


Woy 


cm 


0 . 002 


0.006 


1.5 


1.5 


1.5 


0.8 




cm 


0.003 


. 0018 


0.4 


0.3 


0.6/ 1.2 


0.4 


s 


cm -rad 


0 . 00015 


0.0003 


0.37 


0.07 


0 . 07 


0 . 16 


F 


- 


1 


1 


0 . 026 


0.024 


0 . 24 


0.28 



Table IV-1 Table of High Gain FEL Parameters 
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The following chapter is a description and analysis of 
each PEL in order to characterize them with our simple 
model . 

A. DESCRIPTION OF LCDS AT SLAC 

The Linac Coherent Light Source (LCLS) proposal is to 
use the last kilometer of Stanford Linear Accelerator Center 
(SLAC) to build a PEL able to achieve a 1.5 A x-ray 
wavelength. An electron gun will inject an electron bunch 
that will be accelerated along the last kilometer of the 
SLAC accelerator to 14.35 GeV energy with a peak current of 

3400 A. The normalized emittance is Sj, of 1.5 mm-rad. After 

the undulator, the electron beam will be dumped. The 
accelerator is tunable from 14.3 to 4.5 GeV which 

O 

corresponds to x-ray wavelengths of 1.5-15 A. The electron 
beam particle density is p=2.5xl0^® cm'^ in a beam radius of 
rg=30ja. The undulator will be a linearly polarized with a 

length of L=112 m, period Xo=3 . 0 cm, number of periods 
N=3328, peak magnetic field strength B=1.32T for a resultant 
undulator parameter of K=2.6, and a dimensionless current of 
j=27, 000 . 

The PEL will be a single-pass amplifier starting from 
noise or spontaneous emission; a process called Stimulated 
Amplification from Spontaneous Emission (SASE) . Due to the 
lack of mirrors at this wavelength , numerous proposals are 
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being investigated to guide the light once it leaves the 
FEL. 

B. TESLA FEL AT DESY 

The Deutsches Elektronen-Synchrotron (DESY) in Hamburg, 
Germany, is also planning to develop an x-ray FEL using the 
TeV Energy Superconducting Linear Accelerator (TESLA) . The 
electron beam will be supplied by two 16 km long 
superconducting linear accelerators which were initially 
designed to collide electrons and positrons. 

The specifications for the TELSA FEL are still under 
consideration and the project is being planned in 



conjunction 


with 


the 


LCLS 


with a later 


timeline . 


Many 


parameters 


will 


be 


based 


on experiments 


conducted 


in 


preparation 


for 


LCLS. 


The 


projections are 


based on 


beam 



energy of 25GeV and a peak current of 5kA. 

C. ANALYSIS OF SLAG AND TESLA PARAMETERS 
1. Wavelength 

The purpose of these two lasers is to create a source 
of high power, coherent, hard x-rays at a 1 A wavelength for 
numerous experimental uses. The lack of normal incidence 
optics at this wavelength creates many unique problems. Only 
grazing- incidence optics are possible, and they are still 
under development. The limited optics precludes the use of 
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an oscillator configuration and requires the FEL to operate 
as a single pass amplifier. 

2 . Reducing j 

The dimensionless current j for both x-ray lasers is a 
large number. There are several factors, which lower the 
effective j for actual use in our simulations. 

a) Filling Factor 

In a typical FEL, the radius of the electron beam 
is smaller than that the optical mode. The filling factor is 
the cross sectional area of the electron beam divided by the 
cross sectional area of the optical beam. For LCLS and TESLA 
the cross sectional area of the electron beam is comparable 

to the cross sectional area of the optical beam, so F«1 . For 
LCLS, the radius of the electron beam r^ is about 0.03 mm= 
30 ) 0 .m so the normalized pulse width cr^»0.45. The radius of 

the optical field is approximately W(,»0.02 m= 2 )a.m. 




nini 

Figure IV- 1 Gaussian X-ray Beam and Top-Hat Electron Beam 
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b) Betatron Oscillation 



The number of electron betatron oscillations along 
the undulator in the x-ray FELs is small. The number of 
betatron oscillations along the undulator is inversely 
proportional to the energy of the electron beam. If there 
were a large number of betatron oscillations, then the 
average amount of time electron beam spend within the. 
optical field is of little concern. Betatron oscillations 
cause the electrons to pass from the outside of the optical 
mode, where the fields are weak to the center of the mode, 
where the fields are strong. In order to account for the 
reduced bunching drive of the electrons due to betatron 
oscillations, due to the dimensionless current j is 
decreased by half. 

3 . Rayleigh Length 

The waist size of the x-ray beam is very small with a 
W()=20 jj.m radius. A small Rayleigh waist means a small 

Rayleigh length and large diffraction pattern a large 
distance away. The x-ray beam diffracts away from the 
electron beam contributing to reduced coupling. The Rayleigh 

length, Zq = uwg^ / X , for the LCLS this is approximately 

30m, about one- third the length of the undulator. In a high- 
gain FEL such as LCLS the gain is exponential along the 
undulator. For LCLS, the gain in a Rayleigh length is about 
100m as shown in Figure IV-2. 
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Figure IV- 2 Discrete comparison of gain versus diffraction 

To illustrate the reduced role of diffraction of the x- 
ray beam, divide the undulator into three discrete pieces. 
At the beginning of the undulator the optical power is taken 
to be 1 as shown in Figure IV-2 on the left. In one Raleigh 
length, Zq, along the length of the undulator, the gain is 
100, but only within the region of the electron beam. The 
power in the center is 100, while the x-ray beam outside the 
electron beam has had no gain and remains power 1. The 
strength of the diffracted x-ray beam is only 1% of the 
power of the x-ray beam in the center. In the second 
Rayleigh length, the x-ray within the electron beam region 
has been amplified another factor of 100. The power in the 
center x-ray beam after two Raleigh lengths is 100^. The 
power of the diffracted beam was amplified to 100 in the 
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first Rayleigh length and remains 100. Again this only 1% of 
amplified x-ray beam within the electron beam region. In the 
third length, the power at the end of the undulator has 
become 100^, while the diffracted light has only power of 
100^. The power of the diffracted region is again only 1% of 
the power of the beam in the center at the end of the 
undulator . 

In reality this is a continuous process along the 
undulator. But, just as for this discrete example, the 
continuous gain within the central area of the electron beam 
far out-weighs the gain at the edges of the electron beam. 
To account for this reduction in coupling, the value of j 
can be reduced by half, so that considering betatron motion 

and diffraction the final value is j«6700. This is 
comparable to the other long wavelength FELs in Table IV- 1. 

4 . Emittance 

The emittance for the high gain FEL can be extremely 
important. The x-ray field growth is exponential with an e- 

folding time of 2 /(j / *0.15 for j=5000. For LSCS 

the emittance of ae=0.7, leads to a e-folding decay time of 
w 1 / Gg w 1.4 . Since the decay time caused by emittance 
is greater than the growth e- folding time driven by the 
electron beam, x^ « 1.4 > x^ » 0.1. The net growth rate of 
the x-ray beam is significantly reduced along the undulator. 
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5. Slip 

The electron pulse length is about 186 time longer than 
the slip distance N?t. The extremely short x-ray wavelength 

makes the slippage distance very small and short pulse 
effects are non-existent. Each part of the x-ray pulse 
experiences gain, which is proportional to the local density 
of electrons. 

6. Initial Optical Field 

Since no other laser can radiate near x-ray 

wavelengths, the x-ray PEL must seed itself. The initial x- 
ray field a^ for an x-ray laser is necessarily zero, and the 
laser depends on the SASE (stimulated amplification of 
spontaneous emission) to begin. A simple simulation model is 
to start with a very small initial optical field that is 
representative of the spontaneous emission created in the 
beginning section of the undulator. In the first portion of 
the undulator, the model assumes that an initial optical 
field starts from random fluctuations in the electron beam 
density. The x-ray mode in SASE is known to be linear, 
starting from random noise, and growing through mode 
competition to the desired wavelength. Since the field 
growth in the beginning of the undulator is linear, instead 
of the exponential gain seen later along the undulator, a 
slightly longer undulator is required. 
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D. ELECTRON LASER FACILITY (ELF) 

The Electron Laser Facility (ELF) is an offshoot of the 
Experimental Test Accelerator (ETA) located at Lawrence 
Livermore National Laboratory. The ELF was built to utilize 
the electron beam accelerated by ETA. In each year, 1985, 
1986 and 1987 a paper was published updating that years 
experiments at ELF (references 4,5, and 6). The main 
components of the laser remain essentially the same with 
several modifications from year to year. Here, the three 
phases of the ELF are referred to as ELF85, ELF86, ELF87. 
ELF85 had beam energy of 3.3 MeV and ELF86 and ELF87 had a 
beam energy of 3.5 MeV. ELF85 and ELF86 used a 3m long 
undulator which was extended in ELF87 to 4 meters. The 
undulator for all three is composed of specially shaped 
solenoids with a peak magnetic field of 5 kG. The wiggler 

period remained )t(,=9.8 cm for all three. The number of 
periods was N=3 0 and N=4 0 respectively. The peak current 
initially was 500 A for ELF85, then increased to 850 A for 
ELF86, and finally to 1000 A for ELF87 . The wavelengths 

produced were ^=0.76 cm, 0.72 cm, and 0.175 cm. For each 
experiment, the optical field was contained in a 2.9 cm by 
9 . 8 cm rectangular waveguide and the FELs were operated as 
single pass amplifiers. 
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E. ELF85 

ELF85 has a dimensionless current j=3300. The 
dimensionless beam size is 0.045, and the initial optical 
field of ao=12 was produce by a 60kW pulsed magnetron. 

1. Electron Beam to Optical Field Ratio the Filling 
Factor 

A sketch of the cross section of electron beam and 
optical mode is shown in Figure IV-3 and the filling factor 
is about F=0.03. 
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Figure IV-3 ELF 85 optical field (white) , electron beam(gray) 

diagram 

The electron beam is small enough so that electrons see only 
the peak of the optical beam. There are a large number of 
betatron oscillations along the undulator, unlike LCLS, but 
the electrons continue to see the peak optical field. The 
betatron oscillations increase little effect, since they 
remain within the optical field. There is some radiation in 
the harmonics due to the oscillation, but not enough to be 
significant . 
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2 . Slip 

The electron pulse is 20 slippage distances in length 
for ELF85. This is much smaller than for the LCLS, but still 
much larger than 1. Since the slip is larger than 1, short- 
pulse effects can still be neglected. 

3 . Emittance 

The emittance for the beam was quite large s=0.37 cm- 
rad. This large emittance, assuming a matched beam, means a 
ae=22, and therefore a characteristic decay time of 

Ta=l/o'e=0 . 045 . Gain degradation due to the emittance is quite 
significant . 

F. ELF 86 

The major difference between the ELF85 and ELF86 is the 
modification in the Experimental Test Accelerator (ETA) to 
use a field emission cathode. This increases the peak 
current and makes the pulse longer and the emittance 
smaller. The undulator dimensions remained the same. 



49 



1 . 



Electron Beam to Optical Field Ratio, the Filling 
Factor 
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Figure IV- 4 ELF 86 optical field, electron beam diagram 

The ratio of the electron beam area to the optical beam 
area decreased slightly giving a filling factor F=0.024, as 
seen in Figure IV-4. The electron pulse length remained at 
20 slippage distances with minimal short pulse effects. The 

dimensionless beam size remained nearly the same due to a 
slightly smaller optical wavelength 0.72 cm. The 
dimensionless current increased to j=5000. The input optical 
field is provided by a lOOkW, 35GHz magnetron giving ao=14 . 

The emittance was significantly reduced to s=0 . 07 cm- 

rad. The overall increase in current and better beam quality 
has increase the power output from 80 MW for ELF85 to 180MW 

for ELF86 and the efficiency increased from ri = 5% to r|=6%. 

G. ELF87 

A major change in ELF87 was the use of a high-current 
induction linac. This linac provided a beam energy of 3.5 
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MeV and a peak current of 1 kA. The dimensionless current j= 
1000 and the undulator was also extended to L=4 m. 

Another major change was the shape of the electron 
pulse. The beam length cross section was oval as seen in 
Figure IV-5. However, this beam shape has little effect on 
the model . 

The shift in the optical wavelength to X=0 . 175. cm 
causes the optical mode to become much smaller. The ratio of 
electron beam area to optical mode area is much smaller. A 
process called optical guiding also causes the optical mode 
to remain small. 



a 



U 




Figure IV-5ELF 87 optical field, electron beam diagram 

The net effect is a larger filling factor, F=0.24, and more 
interaction between the light and the electron pulse. 

H. MIT FEL 

The MIT FEL is operated by the Department of Physics 
and Research Laboratory of Electronics, Massachusetts 
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Institute of Technology. This FEL was designed to operate at 
millimeter and sub-millimeter wavelengths. The MIT FEL has 
not been studied as closely, since it has lower gain than 
ELF. The most significant difference between this laser and 
the ELF laser is the undulator. The MIT laser uses a helical 
undulator . 

1. Natural Focusing 

An advantage of the helical undulator is the natural 
focusing of the electron beam. The magnetic field is always 
pushing the electrons toward the center of the beam. Far an 
x-ray FEL, any added benefit from natural focusing could 
result in a significant increase in gain. 

2 . SASE 

The other significant feature of the MIT laser is the 
SASE startup at long wavelengths. From its operation, it can 
be seen that SASE occurs within about the first 20 undulator 
periods. After the first 20 periods, the FEL field growth 
becomes exponential and can be simulated using a small 
initial optical field a.^. 

3 . Electron Beam 

A multielectron gun with a Pulserad IlOA accelerator 
creates the electron beam. The beam is of relatively low 
kinetic energy, 2.3 MeV, with a peak current of 930 amps. 
The pulse length is 7200 times the slippage distance, so 



52 



that the light pulse interacts with the electrons along the 
entire length of the undulator. 

I . CONCLUSION 

In careful study of these FEL's, one can discern which 
parameters have the most influence on the growth rate and 
FEL interaction. A simulation model can often be simplified 
by concentrating on those parameters which create the large 
growth rate or decays, and ignoring those which have no 
discernible effect on the electron/optical field 
interaction. The parameters developed in this chapter can be 
used in a simple model. This is left for further study. 
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V. 



DARMSTADT FEL 



The next section is an analysis of the Darmstadt FEL 
utilizing computer simulations. The Darmstadt FEL is a 

medium gain laser and operates between X=6.6 and 7.8jJ.m 
wavelength, with electron beam energies between 32 and 29 
MeV. It is driven by the S-DALINAC, a re-circulating 
electron linac using superconducting acceleration cavities, 
which ’allow for a continuous wave electron beam. 

A. DARMSTADT PARAMETERS 

The parameters for the Darmstadt FEL are shown in Table 



V-1. The primary 


use 


for the 


Darmstadt 


FEL is to provide 


infrared light 


for 


medical 


research , 


Most of- the 


experiments involve 


medical 


uses such 


as ablation of 


biological tissue. 











Electron Energy-KE 
Peak Current -I 
Electron Pulse Length-lg 
Number of Periods -N 

Undulator period-A-g 
Undulator Parameter-K 

Dimensionless current-j 

Energy Spread-AE/E 
Q-factor 



31 MeV 
2 . 7 A 

600)am 

80 

32mm 

0.78 

1.6 

0.0023 

110 



Table V-1 Darmstadt parameters 



B. BEAM SIZE EFFECTS ON GAIN 

Initially the parameter values obtained for the 
Darmstadt FEL were not properly understood. The electron 
beam radius was given at the end of the undulator, and it 
was mistaken to be the radius at the center of the 
undulator. After several communications and some analysis, 
the discrepancy was discovered. As a result of angular 
spreading along the undulator the beam radius was 0.2mm at 
the waist, and a little less than 1mm at the end of the 
undulator. For an electron beam with Darmstadt's parameters 



and an electron beam 


waist 


radius r^ 


= lmm. 


the 


angular 


distribution would be 


ae=7.3 . 


In the 


time 


it 


took to 



discover the error, an analysis of the influence of angular 
spread on gain was conducted. Despite not having any 
relevance to the Darmstadt FEL, the analysis was interesting 
and may have future relevance. Therefore the results are 
included here. The analysis below was conducted assuming an 
electron beam waist radius of 1mm. 

1. Theoretical Background 

An angular distribution of Gq=7.3 is significant enough 
to degrade the gain of an FEL. A large electron beam radius 
means that there are a significant number of the electrons 
which are injected into the undulator that are off-set from 
the axis of propagation and at an angle with respect to that 
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axis. The larger the diameter, the more electrons that are 
injected off center and at larger angles. The electrons on 
the outer portion of the beam experience a larger magnetic 
field. These outer electrons are "wiggled" harder and 
therefore travel on a slightly longer path. As a result of 
the longer path, the electron slows down in the longitudinal 
direction with respect to an electron in the center of the 
beam. The spread in longitudinal velocities creates a spread 
in phase velocity. A spread in phase velocity then results 
in a reduction in gain. 

The same argument can be used for electrons that are 
injected at an angle with the axis. An electron injected at 
an angle moves toward (or away from) the more intense 

magnetic field off axis. The electron is pushed back to the 

center harder than the other electrons and travels a longer 
distance. The longer distance means a greater spread in 

phase velocity and less gain. 

2 . Method o£ Analysis 

In order to determine the effect of beam size on the 
gain, the effect of value of Ge on the gain is analyzed. The 
gain for each Ge must be calculated at the value of 

desynchronism giving maximum gain. For each Ge, there is a 

different desynchronism value d for which gain is a maximum. 
The amount of time to run enough simulations to determine 

the maximum gains for each Ge would have been prohibitive. 
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Gain versus desynchronism curves were plotted for ae=7.3 and 



ae=3, as shown in Figure V-1. In order to use this 

information to determine how the maximum d depends on ae, 
linear interpolation was used. Assume a linear relationship 
between the maximum desynchronism d and CTe using the form 
d = mOg + b for a straight line. The maximum gain for each 
value of ae can be estimated in subsequent simulations. 

Maximum gain was found for gq = 7 . 3 at d=0.003, and for ae=3 
at d=0.0025. After a small amount of algebra we find that 
m=-0. 00058 and b=0.0066. To test the validity of this 

estimation, ae=5 was estimated to be a maximum at d=0.0037. 

A series of simulations shows that in fact, the maximum 
occurred at 0.0036; in good agreement with the interpolated 
formula. In Figure V-1 it can be seen that there is a flat 
plateau about the maximum of the curve, so our approximation 
for the maximum is accurate enough. In 

Figure V-2 the gain found for each ae, at the maximum d 
determined by the interpolation formula above, is plotted. A 
value of ct 6=0 means perfect injection, no deviation in 
either angle or position. The actual beam radius rg=0.2mm 
corresponds to ae=0.3 which is nearly perfect. The gain for 
ae=0.5 is about 3% less than for perfect injection which is 
an insignificant reduction. 
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Figure V-1 Gain Versus De synchronism 

On the other hand, the original value of ae=7.3 results in a 
77% decrease in gain which is not negligible. In fact, it 
would mean the Darmstadt FEL would not work. 

From Figure V-2 the effects of an angular or radial 
spread on the reduction of Darmstadt gain can be straight- 
forwardly determined. 
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Figure V-2 Gq versus Gain 

The final power P versus desynchronism d curves for <70=7.3 

and 00=3 were plotted, as seen in Figure V-3. The final 
power in strong fields at saturation peaks at small values 
of d where gain is small, and extends out to larger values 
of d . 
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d 

Figure V-3 Power Versus De synchronism 
C . TAPER EFFECTS 

Using the actual Darmstadt FEL electron beam radius of 
rg=0.2mm, and the parameters from Table V-1, the power 
versus desynchronism graphs for several tapers were 
determined. The following positive and negative tapers were 

analyzed: 6=-5.47r, -2.57t, -0.637t, 0, 1.37:, 2.57t, and 5 . 4 t: . 

1. Compiled Data 

Figure V-4 and Figure V-5 are 3-D graphs that depict 
power (on the vertical axis) as a function of both taper 5 
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and desyr.chronis- d. The two figures contain the same da* 
but oresented from different ooints of view. 




Figure V-4 3-D Plot of Power Versus, d, and 5 (Positive View) 

The power versus d curves ?(d) for each value of taper 
6 have similar structures. When d=0, there is no power 

output at any S. At approximately d=0.003, the output power 

is m.aximum for each taper 6. The desynchronisms in the range 
0<d<C.00S do not produce any power for any taper. 
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Figure V-5 3-D Plot of Power, d, and 6 (Negative View) 

The sharp rise in power and high peak power is typical of 



many 


FELs cperaning 


in this 


param.eter regime 


. The 


largest 


peak 


power occurs at 


5=0. 










A high quality 


factor Q 


means there are 


small 


optical 



loses per pass through the undulator. These losses include 
mirror loss and mirror transparency where optical power is 
removed for the end user. A high Q leads to sharper features 
in the desynchronism, curves and to other important features 
explained in the next section. 
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2 . 



Untapered Power Curve Analysis 



To show more detail of the data in Figure V-4 and 
Figure V-5, Figure V-6 and Figure V-9 show power versus 

desynchronism with no taper 5=0 and the highest value of 
negative taper 5=-5.4 ti. 




Figure V-6 P Versus Desynchronism 5=0 
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In Figure V-6, the FEL exhibits limit-cycle behavior 
for de synchronisms between d=0.0008 and d=0 . Oil (ref (2)), 
and thus is not stable. Since the power levels are not 
stable the power in this region is an average of the power 
oscillations . 

3 . Simulation 

Figure V-7 is the output from a single simulation run 
for the Darmstadt FEL. The input variables for this run are 
displayed in the top square. There are nine output windows 
displayed below the input window. 

The top-left window shows the dimensionless optical 

field amplitude |a(z) | at the end of n=6000 passes. The 

middle-left window shows the evolution of the field 
amplitude after at each pass |a(z,n) |, in this case for 
n=6000 passes. The lighter the color, the higher the 
amplitude, as is the case for each middle window. In the 

bottom-left window, the parabolic electron pulse j (z-x) is 

plotted for reference at the beginning and end of the 

undulator; their separation shows the slippage along the 
undulator . 

The bottom-center window plots for reference the weak- 
field gain spectrum as a function of phase velocity, v. The 

middle-center window shows the optical power evolution 
P(v,n) as a function of phase velocity v at each pass n. The 
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cop-center window show the final optical power spectrum P(v) 
at the end of n=600C passes. 

The bottom-right window shows the total optical power 
evolution ?(n) as a function of the pass number n. The 
center-right v/indcw shows the evolution of the electron 

distribution in phase velocity f(v,n) as a function of 
phase velocity and number of passes n. The upper-right 
v.'indow displays the final electron phase velocity 

distribution f(v,6000) at the end of n=6000 passes. 

, ' ' I 

j **** pel Pulse Evolution **** > 

I 

I j=1.6 (7^=1. 1 <70=0 Oq=2.3 6=0.006 | 

Q=110 A=0 A=0 5=0 N=80 5^=0.0001 | 

|a(z,n)| oi^EZSsO.Oe P(v,n) f(v,n) 

Tl 



-4 z 4 -50 V 50 -50 v 50 




-4 z 4 -50 V 50 0 n 6000 

Figure V-7 Trapped Particle Instability 
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4. 



Limit -Cycle Behavior in the Untapered case 



In the untapered undulator, the trapped-particle 
instability occurs for small values of desynchronism at 
d=0.012. In strong optical fields, electrons become trapped 
in the deep potential wells in phase space. These trapped 
electrons oscillate and radiate into side- bands off the 
main wavelength of the optical field. In the center window 
of Figure V-7, the side-bands off the main peak can be seen. 
Due to 'the power in the side-bands, limit-cycle behavior is 
present in the total power P (n) , and can be seen in the 
output power shown in the lower right hand window. The 
output power oscillates between two power levels as the FEL 
oscillates between two frequencies. The trapped particle 
instability does not always cause limit cycle behavior. In 
Figure V-8, the power exhibits chaotic output. 
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**** FEL Pulse Evolution **** 
j=l . 6 ^ ^ 

Q=110 A=0 A=0 5=0 N=80 5^=0.0001 

|a(z,n)l 0liKE^82.5 P(v,n) f(v,n) 




— 4 z 4 -50 V 50 —50 v 50 




-4 z 4-50 V 50 0 n 10000 



Figure V-8 Ctaotic Power Output 

5. Limit-Cycle Behavior in the Negative Tapered Case 

Figure V-9 is the desynchronism curve F(d) for the case 
of negative taper, 6=-5.4r. For this taper 5=-5.47t, limit 
cycle behavior never occurred, or was so small that it was 
not noticeable . 
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Figure V-9 P versus d 6=-5.4 

In Figure V-10, at the specific de synchronism d=0.005 
and negative taper 6=-5.4 we see there are no trapped 
particles. Above a certain value for d, d=0.004, the trapped 
particle instabilities do not occur. This results in smooth 
power output with no side bands around the optical 
wavelength as seen in the center window. The overall effect 
is a lower power output overall but provides for cleaner 
output power. This may be important for some applications 
where a narrow spectrum is important. For other applications 
more power may be important . 
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**** FEL Pulse Evolution **** 
j=l . 6 ^ ^ 

Q=110 A=0 A=0 §=-5.4tcN=80 6^=0.0001 





—4 z 



4 -50 V 50 0 



n 10000 



Figure V-10 No Trapped Particles 

In Figure V-11, the effects of the trapped particle 

instability can be seen in the left and center windows . 
There is no limit cycle behavior as is seen from the steady 
power ?(n) window in the bottom right-hand window. In every 
taper studied, there were values of d for which the trapped 
particle instability occurred but no limit cycles. 
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**** FEL Pulse Evolution **** 



j=1.6 

Q=110 



a^=i.i Gq=0 
A=0 A=0 



! 

Oq-2.3 d=0.003 
5=-5.47CN=80 5^=0.0001 




|a( 



10000 




-4 2 4 -50 V 50 0 n 10000 



Figure V-11 Trapped Particles with Stahle Power 



D . CONCLUSION 

The data collected from simulations of the Darmstadt 
FSL describe how decreasing electron beam quality would 
decrease that observed gain. Also simulations showed that 
there are regions of operation where limit -cycle behavior 
and chaos may be expected to occur in the experiment. If was 
also observed that the trapped-particle instability m.ay 
occur without the presence of limit-cycles. 
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